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Abstract 



Delta degrees of freedom are introduced into antisymmetrized molecular dy- 
namics (AMD). This is done by increasing the number of basic states in the 
AMD wave function, introducing a Skyrme-type delta-nucleon potential, and 
including A^A^ A^A reactions in the collision description. As a test of the 
delta dynamics, the extended AMD is applied to (p,p') recations at E'lab = 800 
MeV for a ^■^C target. It is found that the ratio and the absolute values for 
delta peak and quasielastic peak (QEP) in the ^^C(p,p') reaction are repro- 
duced for angles Giab ~ 40°, pointing to a correct treatment of the delta dy- 
namics in the extended AMD. For forward angles the QEP is overestimated, 
but in generally the agreement between AMD calculations and experimental 
data is good. The results of the AMD calculations are compared to one-step 
Monte Carlo (OSMC) calculations and a detailed analysis of multi-step and 
delta potential effects is given. Along this analysis a decompostion of the 
cross section into various reaction channels is presented and the reaction dy- 
namics is discussed in detail. As important side results we present a way to 
apply a Gallilei invariant theory for (N,N') reactions up to -Eiab ^ 800 MeV 
which ensures approximate Lorentz invariance and we discuss how to fix the 
width parameter v of the single particle momentum distribution for outgoing 
nucleons in the AMD calculation. 
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I. INTRODUCTION 



In the interesting and expanding field of lieavy ion pliysics several models for the treat- 
ment of the reaction dynamics have been developed. Especially microscopic transport models 
which take nuclear mean field effects and two particle collisions into account have been very 
successful in explaining various kind of experimental data. These models are extentions 
of the basic intra nuclear model (INC) [0,0 in which two nucleon collisions are taken into 
account, but the nuclear mean field is not treated self consistently. 

Two major microscopic approaches are used: transport models of the Boltzmann- 
Uehling-Uhlenbeck type (BUU or VUU) for the nuclear phase-space distribution function, 
solving either a non relativistic transport equation (non relativistic BUU) @J^,§] or a rela- 
tivistic transport equation (RBUU) [^] and the quantum molecular dynamics (QMD) []^ or 
its relativistic version (RQMD) p], in which the nuclear phase-space distribution function 
is modeled by the sum of single-nucleon Gaussian wave packets, whose peak positions are 
propagated in time according to the many-body Newton equation, taking into account the 
self consistent mean field as the sum of all two-body potentials. 

In this paper we will discuss especially the antisymmetrized molecular dynamics (AMD) 
0] which is an extension of the QMD and is similar to the fermionic molecular dynamics 
(FMD) proposed by Feldmeier \W\. In the AMD the nuclear wave function is assumed 



to be a Slater determinant of single-nucleon Gaussian wave packets. Other than in the 
FMD the Gaussian width parameter is time independent and, in addition to the pure mean 
field dynamics which is treated in the FMD, the two nucleon collision process has been 
incorporated in the AMD [^j. 

The AMD has proved to be very successful in describing heavy ion collisions at medium 
energies ( <200 MeV/nucleon) [|,|r^-[l^. Because the AMD describes the total system with 
a Slater determinant of nucleon wave packets it has quantum mechanical character, which 
has been demonstrated in the ability of treating shell effects in the dynamical formation of 
fragments. Furthermore it has been shown that ground state wave functions of colliding 
nuclei given by the AMD are realistic and reproduce many spectroscopic data very well 

EM- 

Since the AMD is successful for energies below 200 MeV/nucleon we want to apply it also 
at higher energies Besides the question of Lorentz invariance, which we will also comment 
on in this paper, the opening of inelastic channels, i.e. the excitation of delta degrees of 
freedom, in the elementary nucleon-nucleon collision process is the most important effect 
for i^iab/^ ^ 600 MeV. Therefore the first part of this paper is concerned with the question 
how to incorporate these inelastic channels into the AMD framework. 

Before applying the extended AMD to heavy ion reactions we want to get a good un- 
derstanding of the basic dynamics in our model. We choose to investigate (p,p') reactions 
at delta resonance energies. This is in the spirit of thr research done by Engel et al. pi 



who studied the pion-nucleus reaction to get a better understanding of the pion and delta 
dynamics in the BUU model. Another reason to apply the extended AMD to (p,p') reac- 
tions is that we found in a previous paper that the AMD can reproduce data for (p,p') 
reactions for energies up to 200 MeV very well, so we want to see if this is also true for 
higher energies. 

since 



It is to be noted that the (p,p') reactions have been a major field of interest p3 
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other than in the (p,n) reaction the delta peak is not shifted. Therefore we will investigate 
the delta potential dependence and the reaction dynamics in detail. 

Finally to get a better understanding of the rather elaborated AMD calculation for the 
(p,p') reaction and the experimental data we compare to a one-step Monte Carlo (OSMC) 
simulation for this process. Thereby we can study the influence and importance of mean 
field effects, multi-step processes, and Lorentz invariant kinematics. 

The organization of this paper is as follows: In Sec. ^ we explain the AMD framework, 
discussing the introduction of delta degrees of freedom in the mean field and in the collision 
part; i.e. the delta potential and NN —>■ NA cross sections are introduced. In this section 
we also discuss the first calculations of delta separation energies of 'delta nuclei' in the AMD. 
In Sec. m we give a description of the OSMC which we use for comparison with the AMD 
calculations. As an important basis for the studies presented in this paper, we study the 
question of Lorentz invariance for the (p,p') reaction within the OSMC. We find and show the 
interesting result, that a combination of Lorentz transformation from the laboratory frame 
into the total equal-velocity frame, Gallilei invariant treatment in the equal-velocity frame, 
and following Lorentz back transformation into the laboratory frame ensures approximate 



Lorentz invariant results for energies around i^iab ~ 800 MeV. In Sec. [I^ we discuss the 
results of AMD and OSMC in several steps and compare them with one another and with 
experimental data. First we present the results of the OSMC calculations and compare 
them with experimental data. Before we discuss the improvements which are achieved by 
the full AMD calculation, we explain how to fix all the parameters, especially the width 
parameter u in the single particle momentum distribution of the outgoing nucleons, in the 
AMD cross section calculation independent of the experiment. This is an important point of 
the calculation and should be understood well. After a detailed discussion of the multi-step 
contributions to the cross section in the AMD calculation we proceed to discuss multi-step 
and potential effects by comparing OSMC and AMD calculation with one another. In a 
final step we discuss dependence on delta potential, on elastic nucleon-nucleon cross section 
and on target momentum distribution of the calculations, and possible improvements of the 
AMD calculations. In Sec. ^ we give the summary and an outlook on what should or what 
could be done in future. 



II. FORMULATION OF EXTENDED AMD 
A. Basic idea 

For intermediate energy collisions Eiab/A ^ 600 MeV inelastic channels open in the 
elementary nucleon-nucleon process. The first channel to open is the pion production. As it 
is kown from one-boson exchange model calculations for the elementary process A^A^ — > NNtt 
P^ , |25| the pion production proceeds mainly through excitations of delta states as it is 
depicted in the Feynman diagram in Fig. Diagrams with an excitation of an off-shell 
nucleon are only important for energies near the pion production threshold. 

To include additional processes in the collision part of microscopic models the cross 
sections for this processes are needed. In the pion production case these cross sections 
are well known from experiment, but still there are in principle two ways to include pion 
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production into microscopic models. One way is to treat the A^A^ NNir reaction directly, 
but this approach has several deficiencies. First of all the pions in heavy ion collisions are 
produced in the nuclear medium and therefore can be reabsorbed during the collision. Since 
the pion absorption process in the nucleus is not well understood, it is not clear how to treat 



it after a pion productuion in the nuclear medium. Salcedo et al. calculated density- 
dependent pion absorption probabilities in the nucleus, but these calculations are based on 
local Thomas Fermi approximation and therefore are only valid for ground state nuclei and 
cannot be used in heavy ion collisions. Another drawback of the pure pion and nucleon 
cascade approach is that it is well known from experiment that also delta degrees of freedom 
survive in the nuclear medium. 

Due to this drawbacks we choose another way which is also adopted in other microscopic 
models (see for example Refs. [p],p|, p7| , p8| ). We cut the A^A^ NNn reaction into two steps: 



1) A^A^ — >■ A^A and 2) A — A^vr reaction. This is possible since the decay of deltas into 
pions occupies nearly 100% of the delta width. 

To simulate the produced pion and nucleon momentum spectra correctly we have to take 
a mass distribution of the intermediate excited delta resonance into account. The origin of 
this mass distribution is easily understood by calculating the contribution of the Feynman 
diagram in Fig. to the pion production cross section. In the production cross section a 
term proportional to the square of the delta propagator 

Ga ~ fl) 

^ s - (Ma - ir(s)/2)2 ^ ^ 

will appear, where s is the square of the four momentum of the delta, the delta rest 



mass and r(s) the width of the delta (see Eq. (15)). Therefore off-shell states will contribute 



to the cross section with the weight of a s-dependent factor proportinal to G\. To simulate 
the off-shell character of the delta in microscopic models we populate delta states with a 
mass distribution proportional to G\, interpreting s as mass of the delta, but propagating 
the delta states as on-shell particles (see Sec. |11D| ). 



Now by including also the reverse processes A^A NN and vrA^ — >■ A we get a consistent 
picture of the delta and pion dynamics in heavy ion collisions. Unklike in the pure pion case 
discussed above, the pion absorption process is fixed in this approach, namely, 1) ttN — >■ A 
and 2) A^A — » A^A^ lead to a pion absorption. It has been shown that this treatment of the 
collision term can reproduce many experimental data for the pion-nucleus reactions quite 



well |£T[. Also unlike in the pure pion case the delta is included explicitly as degree of 
freedom. 

As a first step of extending the AMD we discuss in this paper only the extension to 
include delta degrees of freedom. We find that for calculations of (p,p') reactions for light 
targets the treatment of only delta degrees of freedom is sufficient. The decay of the deltas 
is calculated in the final time step of the calculation, when the deltas have escaped the 
the nuclear medium. Later we plan to include also pion degrees of freedom in the nuclear 
medium and the full pion dynamics as explained above. 

Besides the inclusion of the delta in the collision part it also has to be included into the 
mean field part of the theory. Normally this is done by treating the delta as heavy nucleon 
and propagating it in the same mean field as the one for nucleons. More refined calculations 
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have been done which use different mean potentials for deltas and nucleons [^. We adopted 
also a delta potential which is different from the nucleon potential. 

We will proceed in the following Sec. [11 B| by explaining how to describe a system of 
deltas and nucleons in AMD, presenting the basic formulas and also the various potentials 
used in the later calculations. After this we describe in Sec. [11 D| the adopted cross sections, 
i.e. the delta related cross sections, and how we treat the collision term in the extended 
AMD. 



B. Mean field with delta 

The details of the AMD formulation have been explained in Ref. [0. We will only 
comment on the modifications due to the introduction of delta degrees of freedom and those 



points which are important for the understanding of the calculation discussed in Sec. |V, 
In AMD, the wave function of A-nucleon system is described by a Slater determinant 

l<f(^)), 



\HZ)) 



1 



det 



(2) 



where x stands for the spin-isospin function and aj represents the spin-isospin label of the 
j-th single particle state, ctj = p t, p J,, n t, or n |. (pZj is the spatial wave function of the 
j-th single-particle state which is a Gaussian wave packet. 



(r|0z,) 



2z/N 3/4 

1 exp 



TT 



(3) 



where the width parameter v is treated as time-independent in the present work. We take 
z/=0.16 fm~^ in the calculation in this paper. Here Zj is the complex vector whose real 
and imaginary parts, and Kj, are the spatial and momentum centers of the packet, 
respectively. 

For the description of delta states we extend the possible number of spin-isospin functions 
X in Eq. (0). This leads to the introduction of 16 additional states. But since we only use 
spin-independent cross sections and spin-independent delta potential we can reduce the 
number of delta states to four, 



(4) 



or in other words we use a delta spin averaged wave function. For deltas we use the same 
spatial wave functions as for nucleons as given in Eq. (^) leaving also the width parameter 
V unchanged. 

In the same way as in the original AMD the time developments of the coordinate pa- 
rameters, Z = {Zj (j = 1,2, . . . , A)}, due to mean field propagation is determined by the 
time-dependent variational principle; 
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{^{Z)\(^h^ - h)\^{Z)) 
6 / dt ^^rii^rryL = 0' (5) 



which leads to the equation of motion for Z, 

d^ d ($(Z)|iJ|$(Z)) 



Q2 (6) 

where a, r = x, z. 

In the case of the extended AMD the many-body Hamiltonian has to be changed. For a 
A-body system with A^^r nucleons and A^a deltas {A = N]\f + N^) the Hamiltonian is given 
by 

H = T + V 

«JV = 1 «A=1 *^ lN<m 

Nn Na Na 
+ ^^VN^^{iN,i^^) + fAA(«A,jA) • 

It is a sum of kinetic energies for deltas and nucleons and the various two-body potentials 
for nucleon-nucleon, nucleon-delta and delta-delta interaction. 

During the dynamical reaction stage of the collision, the total system can be separated 
into several isolated nucleons and fragments. Since the wave functions of the center-of-mass 
motion of these isolated nucleons and fragments are Gaussian wave packets, each of these 
isolated particles carries spurious zero-point energy of its center-of-mass motion. The total 
amount of the spurious energy of center-of-mass motion can be expressed as a function of Z 



1^,12], which we denote as Esp^s{,Z). The actual Hamiltonian we use in the above equation 
of motion (Eq. (|)) is therefore given by ($(Z)|/7|<I>(Z))/($(Z)|$(Z)) - ^sprs(^). 

As discussed in Sec. Ill A| the deltas in the AMD will have a different mass than the free 



on-shell mass. This is reflected in Eq. by using in the expression for the kinetic 
energy. 

Eq. (P) with definitions Eqs. and are the basic formula for the extended AMD on 
behalf of the mean field, but to do actual calculation the form of the potentials have to be 
fixed. 

For the effective two-nucleon force we adopt the Gogny force which has been suc- 



cessfully used in studying heavy ion reactions with AMD |T2|,[T3|. The Gogny force consists 
of finite-range two-body force and density-dependent zero-range repulsive force. This force 
gives a momentum-dependent mean field which reproduces well the observed energy depen- 
dence of the nucleon optical potential up to about 200 MeV but levels off to be zero at 
higher energies whereas the experimental value at i^iab = 800 MeV is repulsive by about 
50 MeV. The nuclear matter EOS given by Gogny force is soft with the incompressibility 
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K = 228 MeV. Corresponding to the choice of Gogny force, the calculational formula of 



the total spurious center-of-mass energy EspT.s{Z) is taken to be the same as Ref. 1T2|. The 
binding energies of ^^C is calculated to be 92.6 MeV while the observed value is 92.2 MeV. 
The calculated r.m.s. radii of ^^C is 2.55 fm which is reasonable. 

For the delta-nucleon potential we use a zero-range Skyrme-type force dependent on 
nucleon density (pn), 

VNA{rN,rA) = S{rN -rA)[toA + t3Ap]i^{rN)] (8) 
with parameters to a, ^3 a and r. This ansatz leads to a total delta-nucleus energy 

Nn Na 

va = mz)\Y.Y.^NA{^N,^A)mz))/{^z)\^z)) (9) 

= y c/r pA(r) (toAPAf(r) + t3APAf(r)^) , (10) 
where pn and pA are the nucleon density and delta density, respectively. If we take the limit 



of infinite homgenous nuclear matter of Eq. (10), we get the same form of the delta potential 



derived from the delta hole model by Ehehalt et al. in Ref. [^. Therefore we use their 
parameters, toA= —700 MeV fm^, t3A= 1750 MeV fm^ and r = 5/3 in the calculation. The 
parameters are chosen to reproduce the mean field for deltas Ua{po) = dVA^Po) /dpA = —30 



MeV for infinte nuclear matter, since this value is known from pion-nucleus scattering 

It is important to note that the nucleon-nucleon potential we use has finite-range which 
leads to a momentum-dependent nucleon-nucleus potential. On the other side the delta- 
nucleon potential adopted here has zero-range, and therefore the delta-nucleus potential is 
momentum-independent . 

Since there is nothing known about the delta-delta potential we choose f a,a to be zero. 
As a remark we want to point out that the calculation presented in this paper are independent 
of the choice of the delta-delta potential because in the calculated nucleon-nucleus reactions 
there is a maximum of one delta excited at a time. For the future calculation of heavy ion 
collisions we propose to use a delta-delta potential of the same type as the delta-nucleon 
potential Eq. (|) since these potentials should be similar except for spin-isopin factors and 
this choice gives a consistent picture of the interaction. 

Numerically the integrations in Eq. ([I0|) to calculate Va(|$(Z) >) and those needed 
to get its derivatives due to Zj are calculated in a similar way as those described in the 



Appendix A of |]rT| for the density-dependent part of the nucleon-nucleon interaction. 

We also include the Coulomb potential in our calculation, but the results are not sensitive 
to this. 



C. Delta nuclei 



To construct the ground states of colliding nuclei in the AMD the frictional cooling 



method [p!7H2(| is used which leads to a consistent picture of mean field dynamics and 
ground state nuclei. It has been checked that wave functions given by AMD are realistic 
and reproduce many spectroscopic data very well. 
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As the first check of the incorporated delta mean field dynamics we calculated the binding 
energies of 'delta nuclei'. This is in a way an academic study since no 'delta nuclei' exist or 
can be detected, but it gives us information on the kind of wave functions we can have in a 
delta-nucleus system. In future this kind of wave functions could be used for studying the 
admixture of delta states in ground state wave function [p2| . 

We calculated the delta separation energies for 'delta nuclei'. For this we used the 
frictional cooling method to cool first an ensemble of protons, neutrons and one additional 
delta. Then we cooled the same ensemble of nucleons without delta. The difference of the 
minimun energies of these two calculations is the given delta separation energy. The results 
for delta separation energies with the parameters given in Sec. [11 B| are shown in Fig. 0. 



The results in Fig. ^ can be understood in the following way: The delta separation 
energy is determined by the overlap between the delta Gaussian wave function and the 
residual many-nucleon system, as can be seen in Eq. ([I0|) . Therefore the minimum energy 
state calculated with the frictional cooling method leads to a delta state with maximum 
overlap with the nucleons. Because of the 3/2 isospin nature of the delta the configuration 
of maximum overlap is Pauli allowed. The mass dependence we get in Fig. ^ for the delta 
separation energy reflects the structure of the ground state wave function in AMD, which 
was able to describe the alpha clustering of light nuclei. We found that the delta separation 
energies for ^He and ^Be are the same. This can be understood on the basis of the two 
alpha clustering of the ^Be ground state. In the ground states of both nuclei the delta is 
in the middle of an alpha cluster, since this is the state of maximum overlap of delta and 
nucleus wave functions. Therefore they have the same delta separation energy. Because of 
the disolving of the alpha structure and the existence of larger zones with higher nuclear 
densities for heavier masses, also the delta separation energies increases for heavier masses. 

D. Deltas in the collision term 

The idea how to incorporate delta resonances in the collision part of microscopic theories 



was already discussed in Sec. |II A| . Here we give the cross sections adopted in the AMD. To 



include the deltas in the collision term we used the same cross sections described in pT| , p8 
therefore we will only recall the most important formulas. 

1. NN NA reaction 

The cross sections for the A^A^ —>■ NA reaction are derived from the VerWest-Arndt 
parametrisation for pion production cross sections in Ref. Asuming that all pions in 



the A^A^ — >■ NNtc reaction are produced via a delta resonance we get 

p + p-^n + A++ (710 + 1/2^11 
p + p^p + A+ 3/2aii 



n + p-^p + AO l/2aii + l/4aio (-.-.. 

n + p^n + A+ 1/2^11 + l/4aio 

n + n^p + A" (Tio + l/2crii 

n + n^n + A^ 3/2 an , 
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for the cross sections where cTj/ denotes the cross sections used in the VerWest-Arndt 
parametrisation for given isospin (/) and (i) of the final and initial nucleons, respectively 
. The angular dependence of the NN — > NA is chosen in the same way as in Ref. . 



As discussed in Sec. |1I A| it is not enough to define the cross section for the NN NA 



reaction but we also need to use an appropriate mass distribution for the delta resonance. 
We choose 

1 r(M)/4 



ttMa (M - Ma)2 + r(M)V4 

(14) 



with the delta rest mass Ma=1232 MeV and the momentum- dependent delta width 

'qV fvjq) 
^Qr) M [viqr) 



m)^(^X^m)'r. (15) 



with 
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/32 + g2 ' 
r^ = 110MeV , 
P = 300 MeV , 

where q and denote the pion momentum in the rest frame of the delta for delta mass M 
and on-shell mass M = Ma, respectively, defined as 

To get information on the validity of the adopted delta mass distribution we calculated 
the neutron momentum spectrum in the p(p,n)p7r"'" reaction and compared to the experimen- 
tal shape of the spectrum. In Fig. ^ the result of this calculation is shown. The experimental 
data is a measurement of the neutron spectrum for laboratory angles 6„ = 0° — 6° in co- 



incidence with a proton and a tt"'" for proton energy £^iab = 830 MeV [^. This reaction 
proceeds mainly via an excitation of a A"*""*" state. The calculations which are shown were 
made with OSMC explained in detail in Sec. |IT| using a proton target: A delta excitation 
in a proton-proton collision is simulated choosing the mass and momentum according to 
the given mass distribution (Eq. (|T^) and the scattering angle according to the differential 
cross section of Ref. , respectively, by Monte Carlo method. The deltas are assumed to 



decay isotropically in their rest frame. Taking into account the proper transformations of 
the momenta we calculate the cross section for neutrons emitted in coincidence with a vr"*" 
and a proton. 
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Besides the calculation with the mass distribution of Eq. ([M]) we made calculations with 
a mass distribution multiplied by a phase space factor 



/phasespace(^^) ^ p^(^M)"^''''^ ^^'^^^ 

with 

pi(M)=l'-<'""--'^']l'-''"" + ^^''' . (18) 

where s is the square of the four momentum of the system. pa{M) is the momentum of 
the delta and pa{M)"^°'^ is the possible maximum momentum of the delta assumed for 
M = nin + m^r. With this phase factor (Eq. ([T7|)) we take into account the fact that the 
final phase space for maximun delta mass, M = ^/s — ruN, is zero. In Fig. ^ calculations 
for the p(p,n)p7r+ reaction with n = 0, n = 2 and n = 3 in Eq. ( jlT]) are shown. All 
calculations are scaled with the same factor to reproduce the experimental data. By taking 
into account only the mass distribution without phase space suppression factor, n = 0, solid 
line, we find that the neutron spectrum is not reproduced so well. Too many heavy A^^ 
with, consequently, low momentum neutron partners are formed in proton-proton collisions. 
Therefore the neutron momentum spectrum is shifted to lower energies. On the other hand 
using a strong suppression factor, n = 3, leads to an underestimation of low momentum 
neutrons as seen in Fig. ^ (dashed line). The best fit for the neutron spectrum in the 
p(p,n)p7r+ reaction with the OSMC simulation is attained when using tt, = 2 in Eq. ([T7|) for 
the suppression factor as seen in Fig. ^ (dashed dotted line). Both low and high neutron 
momenta can be described rather well. In Sec. |IV A| we will comment on the dependence of 



the (p,p') reaction on the choice of the mass distribution. 



2. NA NN reaction 



Since there is no delta beam or delta target available the cross section for the AA^ NN 
process has to be derived by theory. It has been pointed out by several authors PB| , P7|PT[| 
that the naive detailed balance formula is not appropriate for this process. The simple 
argument to understand this point is that in the nucleon-nucleon collisions in which the 
deltas are formed the cross section has to be folded with the mass distribution to get the 
mass of the delta in the simulation, as described in Sec. IIP 1 , but in the NA NN 
reaction the delta in the simulation has a definite mass and therefore the mass distribution 
is not present for the calculation of the A^A —>■ NN cross section. For a more detailed 



discussion of this problem see Ref. |21 



37[] , which was proven to 



For simplicity we adopted the formula derived by Wolf et al. 
give a realistic description of the A^A — »• A^A^ cross section (see discussion in Ref. 
Therefore we use 



4 Pa 



oo 

(mjv- 



F{M^)dM^ 



{y/s-nif^)'^ 

(mjv+m^)2 



F{M^)dM^ 



(19) 
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with the <7nn->na cross sections from VerWest-Arndt and the mass distribution F(M'^) 
given in Eq. (0). In Eq. (0) and are the initial deha and the final nucleon momenta 
in the center of mass frame, respectively, whereas the factor 1/4 is due to spin averaging 
and a symmetry factor for identical particles in the final state. Since this cross section (Eq. 
|T9|) is infinte for zero energy coUisons we apply a low momentum cut off of a"^""^ = 100 mb, 
which takes the screening of the delta in the nuclear medium into account. In case of the 
(p,p') reaction the screening is not effective because the produced deltas are fast deltas. 



3. NN elastic scattering 



For the nucleon-nucleon elastic scattering we used a new parametrisation of Cugnon 
which is isospin dependent and therefore takes into account proton-neutron and proton- 
proton collisions differently |^ . This is not the case in the old isospin-independent Cugnon 
parametrisation [|I|, which is mainly adopted in other microscopic models, in which all 
isospin channels are treated the same way, in contrast to experiment. The isospin-dependent 
parametrisation also includes a back-ward peak for the proton-neutron scattering which is 
not included in the isospin-independent parametrisation. 



E. Simulation 



The actual simulation of the AMD is done in several steps. In the first step the ground 
states of the target and projectile nuclei are calculated by the frictional cooling method. 
With these ground state configurations the collision is calculated in the CM frame (see 
discussion about Lorentz invariant scattering in Sec. [Ill B]) using the above described cross 
sections and mean fields for the particles. The delta decay is not yet incorporated in the 
AMD therefore deltas will survive untill the final state of the calculation. The decay of these 
deltas is calculated asuming an isotropic decay in the rest frame of the delta. In this step the 
relativistic expression for the pion energy is used. In the final step the statistical decay of 
primordial formed fragments is calculated. Primordial fragments mean the fragments which 
are present when the dynamical stage of the reaction has finished. These fragments are not 
in their ground states but are excited, and they decay through evaporation with a long time 
scale. In this paper, the switching time from the dynamical stage to the evaporation stage 
was chosen to be 150 fm/c. Statistical decays of fragments were calculated with the code of 
Ref. P5| which is similar to the code of Piihlhofer |HD|. 



F. Calculation of cross section 



The experimental data we want to compare with in this paper are given in the labo- 
ratory frame. As discussed in Sec. [llllj| we choose a combination of Lorentz and Gallilei 
transformation to achieve approximate Lorentz invariance using a Gallilei invariant theory. 
The combinations of transformations has to be also taken into account for the cross sec- 
tion calculation. In the equal-velocity frame the Gallilei invariant expression of the double 
differential cross section has the following form: 
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dVtdp Jo dVtdp 
(iW(p, h) 



protons 

P^{P) = l(p|0Z,)r = (-^)'^%xp[--^(p-K,)2] (21) 



where d^p = p'^dpdQ and ( )b stands for the average value over the events with impact 
parameter b. In this formula, the outcoming protons are expressed by Gaussian wave packets 
with momentum width h^/u. 

Applying the Lorentz transformation into the laboratory frame we get 

r2nbdb( J: p',{p'))p" (22) 

dn'dp^ Jo \i=i.^ated 

protons 

with 

P\{P') = P^{P)§J (23) 

where p and E are the relativistic expressions for the momentum and the energy in the 
equal-velocity frame, p' and E' the relativistic momentum and energy in the laboratory 
frame and Pi(p) the function given in Eq. (pl|). 



III. ONE-STEP MONTE CARLO SIMULATION FOR (P,P') 

A. Formulation 

In order to get a better understanding of the rather elaborate AMD calculations, we also 
perform one-step Monte Carlo (OSMC) simulations of the (p,p') reaction. In this section 
we explain the OSMC simulation. In the OSMC simulation the summation over the impact 
parameter / 2T[bdb is not made for simplicity. The projectile protons are initialized in the 
laboratory frame according to the relativistic momentum energy relation. Since we studied 
the influence of the target momentum distribution on the results for the (p,p') reaction 
we initialized different target momentum distributions. As one approach we attributed a 
random momentum smaller than the Fermi momentum at normal nuclear density |p| < |kf | 
to the nucleons. As another approach we used the local Thomas Fermi ansatz to simulate 
the momentum distribution by first choosing a position in coordinate space according to the 
space density 



P{r) = Po 



l + a{ -V 

a/ 



exp 



r y 
a J 



(24) 



which is based on the harmonic-oscillator model, with parameters a = 1.247 and a = 
1.649 fm for the ^^C target determined through experiment HTI]. And then we attributed a 
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momentum randomly smaller than the local Fermi momentum, which gave us the momentum 
distribution f{p), 




/W^f *4..VW|||^ .kAr)-[^p(r)) , (25) 

where Q{x) is the step function; 9(x) = 1 for a; > 0, 0(x) = for x < 0. 

As next step we transform the energy and the momentum of the projectile proton from 
the laboratory frame into the equal-velocity frame of target and projectile with Lorentz 
transformation. The transformation f] for incident proton in z-direction and a target with 
A nucleons with four momenta 

p^ = {Ep,0,0,pJ (26) 
p^ = {Am„, 0,0,0} (27) 

(28) 



into the equal-velocity frame with condition 



/ / 4 

^ _P£_ 

^; e'a 



(29) 



IS 



{/3.,/3„/3,} = {0,0,^-./(^)2-l} , (30) 

Pz V Pz 

where A is the number of target nucleons and m„ the nucleon mass. It can be seen that 
this transformation is not dependent on the number of target nucleons. This is natural since 
the total equal- velocity frame is the same as the two nucleon CM frame, if we disregard the 
Fermi motion of the target nucleons. 

We choose to tranform into the total equal-velocity frame, because in this frame non 
relativistic approximations can be applied. In the equal- velocity frame we change to non 
relativistic kinematics using the obtained momenta and the non relativistic expression for the 
energy. Also for the transformation of the target nucleon momenta into the equal-velocity 
frame we use the non relativistic Gallilei transformation. In the next section it will be shown 
that this kinematical treatment of the scattering, which we also adopt in the AMD, leads 
to approximate Lorentz invariance. 

The two nucleon collision in OSMC is calculated in the CM frame of the two nucleons 
with Gallilei invariant kinematics. As possible collision processes we incorporated delta 
excitation and elastic scattering. We fixed the ratio of delta excitation to elastic events 
beforehand, which is well justified, since this ratio is mainly determined by the initial energy 
of the proton. The cross section ratios for i?iab = 800 MeV which are incorporated in the 
calculation are given in Fig. ^ 

In case of a delta event we have to determine the mass of the delta and, since we take 
into account the charge of the particles, the charge of the delta. The charge and mass are 
chosen according to the ratio given by the cross sections and the mass distribution in Sec. 
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II D 1| , respectively. To calculate the cross sections the ^/s of the two nucleon system is 
needed. This we calculate in the two nucleon Gallilei CM frame using the non relativistic 
expression for the energies. In the next section the validity of this evaluation is discussed. 

The three momentum of the outgoing particles, deltas or nucleons, are determined ac- 
cording to the angle-dependent cross sections in the two nucleon CM by Monte Carlo method. 
Since in the delta excitation kinetic energy is transformed into mass energy the final relative 
momentum in the CM frame for the NN — > A^A reaction decreases. 

The decay of the excited deltas is calculted in the delta rest frame assuming an isotropic 
decay. For the delta decay we use a Lorentz invariant prescription because the pions are 
light particles and need to be treated in relativistic invariant kinematics. It is important 
to mention here that the Lorentz transformation in the delta decay is only adopted to 
transform from the delta rest frame to the equal-velocity frame of the total system. The 
charge of the pions and nucleons are chosen according the probabilities given by the Clebsch 
Gordan coefficients for the decay. 

To calculate the multiplicity distribution of the outgoing protons or neutrons in the 
laboratory frame for given momenta and angle, the momenta of the particles are back trans- 
formed from the total equal- velocity frame to the laboratory frame by Lorentz transformation 
with Phack = —P given in Eq. (^). The multiplicity distributions are calculated with the 
formula 

d'^N N- 

(31) 



dpidVti Ap27r A cos QNtp 

where Ni is the number of nucleons with momentum pi — Ap/2 < p < Pi + Ap/2 and cos Oj — 
A cos G/2 < cos Gj < cos Gj + A cos G/2, and Ntp the number of performed simulations. In 
the calculations shown we use Ap = 0.02 GeV and AG = 2.99°. In the calculation of cross 
sections for the (p,p') reaction we use a factor Xfactor to scale to the experimental data: 

-Xfactor ■ (32) 



dpidVti dpidVti 

The factor Xfactor has a dimension, which is because we do not integrate over the impact 
parameter. The number of simulations was varied according to the accuracy needed. 

B. Lorentz invariance 

To justify the approach used for the scattering kinematics in the OSMC simulations and 
also the AMD calculations we compared the OSMC simulation described in Sec. [Ill A| with 



simulations using pure Gallilei and pure Lorentz invariant kinematics, respectively. For this 
comparison we did one-step elastic scattering simulations using an isotropic angular distri- 
bution for the scattering in the nucleon-nucleon CM frame disregarding the charges of the 
particles. The results for the multiplicity distribution for Lorentz invariant kinematics and 
OSMC simulations are shown in Fig. ^ The momentum of the target nucleons is chosen 
randomly between zero and the Fermi momentum kf = 270 MeV. The calculations shown 
are for i^iab = 800 MeV and scattering momenta in the laboratory frame. We show calcula- 
tions for laboratory angles G = 15° and G = 60°. The solid line represents results with the 
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Lorentz invariant kinematics and the dashed hne results with the OSMC. The results agree 
very well with each other. In Fig. ^ the results of Lorentz invariant calculations (solid line) 
are compared with Gallilei invariant scattering calculations. Again the calculations were pre- 
formed for laboratory angles G = 15° and G = 60°. The Gallilei invariant calculations differ 
significantly from the Lorentz invariant results. We did calculations starting with correct 
relativistic energy (dashed line in Fig. and correct relativistic momentum (dashed dotted 
line in Fig. Both calculations deviate strongly from the correct relativistic calculation 
(solid line in Fig. ^). 

To show that the two-nucleon ^^-distribution used in the OSMC simulation and the 
AMD agrees with the relativistic prescription we compared the y/s distributions of the 
OSMC and the Lorentz invariant scattering. Again there are only small differences for 
the resulting y^-distribution as seen in Fig. ^ OSMC simulation (dashed line), Lorentz 
invariant kinematics (solid line). The averaged deviation of ^/s is less than 1% and, since 
the elastic and inelastic nucleon-nucleon cross sections in the considered energy region have 
a smooth dependence on a/s, the difference between OSMC and purely Lorentz invariant 
kinematics is not important. 



IV. COMPARISON WITH EXPERIMENT 



As first step to understand the employed delta dynamics in AMD we have calculated 
the inelastic scattering of protons on light nuclei. We have compared the calculation with 
experimental data in which both the delta excitation and the quasielastic peak (QEP) are 
prominent. In the following discussion we will see that behind the seemingly easy to under- 
stand experimental data complicated and interesting processes are hidden. 



A. One-step Monte Carlo calculations 

First we want to discuss the results of the OSMC calculations. The total reaction prob- 
abilities for the (p,p') reaction for i?iab = 800 MeV incorporated in the OSMC are shown in 
Fig. |. The ratio between inelastic to elastic cross section is approximately 2 to 3, which 
points out that inelastic channels are very important for this energy region. Deviations from 
the total reaction probabilities in the differential cross sections discussed in the following, 
are due to the anisotropic cross sections used in the OSMC and due to kinematical reasons. 

In Fig. ||we show the OSMC calculations (solid line) for the (p,p') reaction at ii^iab = 800 
MeV for laboratory angles G = 15°, 30°, 40° and G = 60°. In the experimental data for 
G = 15° the QEP can be clearly seen and the contribution of inelastic or multi-step scattering 
to the cross section is well separated. In Fig. |] we show two different experimental data 
for G = 15°, which differ only in the QEP region. This difference is claimed to be due 



to the resolution of the experimental apparatus [|^. We scale the OSMC calculation to 
reproduce the momentum region p'=600-800 MeV at G = 15°. To achieve this, we have 
to use Xfactor = 260 mb which is larger than the measured total cross section (209 mb for 



Eiab = 860 MeV in Ref. [||] and 229 mb for E^^ = 553 MeV in Ref. ). 

We find that in the OSMC calculation the QEP for G = 15° is largely overestimated and 
for G = 30° and G = 40° the QEP is much more pronounced than in the experimental data. 
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The high momentum tail in the experimental data at O = 15° is physically not allowed and 
consequently cannot be reproduced. 

For G = 30°, 40° and G = 60° there is an underestimation of the low and high momentum 
tail of the cross section in the shown momentum region. The underestimation is especially 
strong in the momentum region p' = 500-800 MeV for G = 30° and p' = 500-700 MeV for 
G = 40°and . Another defect of the OSMC calculation is that the flat behaviour of the cross 
section for p' = 1000-1200 MeV at G = 15° is not reproduced. Instead there is a dip in the 
calculated cross section. 



B. AMD calculations 

Before discussing the results of the AMD calculation we have to comment on the cross 
section calculation in AMD for the (p,p') reaction. A straightforward approach leads to Eq. 
(P^) for the double differential cross section. This formula is based on the assumption that 
every nucleon in AMD has the same width z/, but this is not true for all nucleons in the 
(p,p') reaction. For example it is obvious that the incoming proton has a definite energy 
and therefore zero momentum width v. Therefore the role of the width parameter has to 
be reconsidered. Concerning the width parameter v in the AMD, we recall, that the width 
parameter u in the AMD is improtant for reproducing the total density and momentum 
distribution of ground state nuclei. In the collision prescription in the present version of the 
AMD (see Ref. [[l^]), only the Gaussian distribution in the coordinate space is taken into 
account, whereas in the momentum space only the centers of Gaussian wave packets are used. 
Therefore in nucleon-nucleus collision calculation with AMD the momentum distribution of 
the target is only partly explored dynamically. For heavy ion reactions the effect of the high 
momentum tail of the wave function is negligible as long as subthreshold processes are not 
important for the calculated process, but, as will be become clear in the following, the QEP 
region of the (p,p') reaction is sensitive to the width parameter z/ of the outgoing nucleons. 

In order to solve the above sketched problem, for getting a reliable exploration with AMD 
of the momentum distribution of the target in the (p,p') reaction, two lines of argumentation 
can be followed. One is to take into account the momentum distribution of the colliding 
nucleons in the collision prescription, and the other is to change the width parameter v of 
the outgoing nucleons. In the first approach the width parameter v of the incoming proton 
should be zero, at least in the simulation of the nucleon-nulceon collison. Technically this 
can be taken into account in the collision prescription, but it is not clear what width should 
be attributed to the incoming proton after the first collision. We follow the second approach 
to get a reliable treatment of the momentum distribution of the target, since there is a way to 
determine the width v of the outgoing nucleons independently of experimental data, namely, 
the AMD calculation should reproduce the width of the QEP caused by the Fermi motion. 
In course of doing so we did calculations both with the OSMC and a modified version of the 
AMD, the 'one-step AMD'. In case of the OSMC we used both, a box momentum distribution 
for the target, like in nuclear matter, and a Thomas Fermi momentum distribution for the 



target, as explained in Sec. [Ill A| . In the 'one-step AMD' calculation we allowed only one 
possible collision in each event after which the particles propagated in their mean field. 
We determine the width v of the outgoing nucleon so that the 'one-step AMD' calculation 
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reproduces the width of the QEP given by the OSMC calculation with the Thomas Fermi 
momentum distribution for the target. 

As can be seen in Fig the QEP at O = 15° is much narrower than at G = 30°. The 
width parameter v in the 'one-step AMD' calculation is fitted so as to reproduce width of 
the QEP caused by the Fermi motion. In the actual calculations we did not change the 
width parameter v but we used a cut of the Gaussian momentum distribution 

= for [^(p - K,)2] > g,,t , (33) 

and renormalized the function p'(p) in Eq. which has the same effect as adjusting the 
width parameter v in the single particle wave functions. 

We found that in the 'one-step AMD' calculation for ^^C(p,p') reaction we have to use 
Qcut = 0.5 for e = 15°, Qcut = 1-0 for 6 = 30°, Qcut = 1-6 for 6 = 40° and no cut for 
= 60° in order to reproduce the width of the QEP in the OSMC simulation with Thomas 
Fermi momentum distribution of the ^^C target. 

In all AMD calculations we use the above determined cuts. 

By comparing the 'one-step AMD' with OSMC calculations with box momentum distri- 
bution, the nuclear matter case, we find that we only have to apply a cut at G = 15° to 
reproduce the effect of the Fermi motion. This is a nice result which states that for cross 
section calculations for heavy ion collisions no cuts are needed, since large amount of nuclear 
matter is involved in the reaction. 

After fixing all parameters we can proceed to discuss the results of the AMD calculations 
which are shown in Fig. ^ (solid line). In general we find a better reproduction of the experi- 
mental data with the AMD than with the OSMC. Here we will only state the improvements 
established by using the AMD instead of the OSMC. In the following section (Sec. |1V Cp 



we study the reaction mechanism in the AMD calculation in detail and in Sec. |IVD| we will 
discuss the origin of the improvements in the AMD. 

Instead of an underestimation of the low momentum region for angles O = 30° and 
O = 40° in the OSMC, we find that the experimental data can be better reproduced in the 
AMD calculation. In general the agreement with the experimental data for these angles is 
much better in the AMD than in OSMC, because the smooth behaviour of the experimental 
data for O = 30°, and especially for B = 40°, is much better reproduced. In the case 
of the OSMC we do not get absolute cross sections and therefore the comparison of the 
calculations with the shape of the momentum (p') dependence is important and this is much 
better reproduced in the AMD calculations. 

The AMD can also reproduce the high momentum tail of the cross section for O = 30° 
better and for O = 40° and for = 60° much better than the OSMC. Also the the flat 
behaviour of the cross section for p' = 1000-1200 MeV at = 15° is reproduced. 

So we find that the agreement between AMD calculation and experiment is much better 
than the one between OSMC and experiment. The experimental data can be quantatively 
reproduced in AMD, whereby the elastic processes seem to be overestimated in forward 



direction. In Sec. P^VE| we will discuss how to improve the AMD results, but in the next 
sections we will first try to get a clearer understanding of the reaction dynamics in the AMD 
for the (p,p') reaction. 
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C. Multi-step contributions to the cross section 



To get a better understanding of the reaction dynamics we calculated the decomposition 
of the cross section into different multi-step contributions and reaction processes. The 
result of this decomposition is shown in Fig. We show calculations for four angles: 
= 15°, 30°, 40° and = 60°. The elastic multi-step contributions shown in Fig. ^ 
are protons with one-step elastic scattering (ELA-lstep, long dashed line) and protons with 
two-setp elastic scattering (ELA-2step, dashed line). Nucleons involved in inelastic collisions 
can contribute to two different sources of outgoing protons: protons emitted as delta decay 
products and protons which acted as collision partners in the delta excitation process. The 
contributions of these different sources in inelastic collisions are also shown in Fig. |l^: 
Protons emitted as partner of an excited delta in the first chance collision (A-lstep(N), 
dotted line), and protons from delta decay after the first chance collision(A-lstep(A), dashed 
dotted line). 

In the calculation for 6 = 15° we see that mainly the one-step scattering process (long 
dashed line) contributes to the QEP. The two step contribution is much smaller for all 
momenta. A similar behaviour is found for B = 30°, where the two step contribution 
(dashed line) gets more important but is smaller than the one-step contribution for almost 
all outgoing proton momenta. At = 30°, G = 40° and B = 60° the two step contribution 
dominate the high momentum tail of the proton spectrum, since one-step scattering cannot 
reach this region of the phase space easily. 

For = 30° and = 40° we find that for the lower momentum part of the cross section, 
below the QEP, many higher step processes contribute to the cross section. Besides the 
processes shown in Fig. ^ two-step and three-step elastic processes followed by an inelastic 
process contribute to the cross section. Therefore the reproduction of the smooth behaviour 
of the cross section in this momentum region is quite sensitive to the ratio of inelastic to 
elastic cross section and the multi-step contributions. 



D. Comparison of AMD and OSMC: Multi-step and potential effects 

After having discussed the reaction dynamics in the AMD in detail we can comment on 
the differences between AMD and OSMC and what kind of physical effects are important 
for the (p,p') reaction. 

As one could have expected, the main differences of OSMC and AMD are due to multi- 
step processes, but the details of the origins of the differences are much more interesting 
than this plain statement, and therefore we discuss this question in detail. Also we have to 
investigate the role of potential effects which could interfere with the multi-step contribution 



effects. To make the above statements more transparent we show in Fig 11 the outcome of 



the 'one-step AMD' calculation (dashed line) we mentioned in Sec. [IV B|, in which we did 



not allow anymore collisions after the first proton- nucleon collision ocurred in the AMD, 
propagating the particles from this time step onwards only in their mean fields without 
collisons. Therefore the only difference between 'one-step AMD' and OSMC is the treatmant 
of the mean fields in the AMD and the differences between 'one-step AMD' and AMD are 
only due to multi-step contributions. The results of the 'one-step AMD' should not be 
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confused with the one-step contributions to the full AMD calculations discussed in the last 
section. 

First we want to discuss the differences due to multi-step contributions and the improve- 
ments achieved by this in reproducing the experimental data with the AMD. Therefore 
we compare 'one-step AMD' (dashed line) and AMD (solid line) calculations in Fig. ITT] on 



the basis of the multi-step decompositon of the AMD results shown in Fig. |10|. Multi-step 
contributions which involve at least one inelastic step are mainly responsible for a better re- 
production of the experimental data for momentum below the QEP at O = 30° and G = 40° 
in the AMD calculation. The high momentum tails above the QEP at (G = 30°,) G = 40° 
and G = 60° are better reproduced in the AMD by mainly multi-step elastic scattering. 

The apparent quantitative agreement of the heights of the QEPs for G = 30° and G = 40° 
in the 'one-step AMD' and the AMD is not a trivial result. The reduction of the QEP due to 
multi-step scattering is balanced out by an enhancement due to multi-step processes in the 
AMD. At G = 60° multi-step processes lead even to an increase of the cross section in the 
QEP region in the AMD calculations. For G = 15° the net effect of multi-step contributions 
is a reduction of the QEP. 

Another important improvement of the AMD calculations due to multi-step contributions 
is that the fiat behaviour of the experimental data at G = 15° for momenta p'= 1000-1200 
MeV can be described unlike in the OSMC where there is a dip in the cross section. 

Concerning the potential effects in AMD we find no significant difference between 'one- 
step AMD' and OSMC calculations. There are no apparent differences caused by the in- 
clusion of potentials for the particles. (See more detailed discussion about delta potential 
effects in the next section.) 



E. Dependence on delta potential, cross section and target momentum distribution, 

and improvements 

Since we have now established a detailed understanding of the reaction dynamics in the 
AMD we can study the dependence of the results on the physical input quantities and what 
we can learn about these in the (p,p') reaction. 

One important question which is always asked, but not yet well understood, is the 
question of the delta potential and the behaviour of the delta in the nuclear medium. We 
want to comment on this from the viewpoint of the AMD calculations made in this paper. 
In the previous section we stated that there is no apparent potential effect when comparing 
OSMC and AMD results, but here we will investigate the potential effects in more detail. 

To study the effect of the delta potential we did calculation with a potential which has 
the same density dependence as the one in Eq. ||, but is shifted to reproduce ^7a(po) = —80 
MeV. This choice of the potential is motivated by the following idea: We assume that the 
delta potential is momentum-independent and has the original value U^{pq) = —30 MeV. 
We recall here that we have not treated the momentum dependence of the nucleon potential 
correctly in the present version of the AMD and our nucleon potential at i^iab = 800 MeV 
is Un{po) ~ MeV. If we used the correct value Un{po) ~ 50 MeV, it would lead to the 
difference Un{po) — Ua{po) ~ 80 MeV. Hence using a delta potential with U/\{po) = —80 
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MeV in the present version of the AMD leads to the desired value Un{po) — Ua{po) ~ 80 
MeV. 

The AMD calculations with shifted delta potential are shown in Fig. ^ (dotted line). 
There is no apparent difference between the calculations with different delta potentials. 
Only if we calculate the multi-step decomposition of the cross section we find that the high 
momentum threshold of the one-step inelastic contributions are different for different delta 
potentials. But these differences cannot be seen in the sum of all contributions because in 
the high momentum threshold region for these processes their contributions to the cross 
sections are just a few mb and therefore are hidden below the multi-step and elastic step 
contributions to the cross section. This is also the reason why we did not find any apparent 
potential effect in the previous section. 

We point out at this point that the choice of the delta mass in the A^A^ NA process is 
not treated fully self consistent in the present version of the AMD. To determine the mass 
of the delta we use the delta mass distribution Eq. |TB| dependent on the -^i of the initial 
nucleon-nucleon pair, but we do not take the effect of the delta potential on the delta mass 
distribution into account. This improvement should be made in future AMD calculations. 
Still the above findings about the influence of the delta potential on AMD calculations for 
experimental data for the ^^C(p,p') reaction would hold. This is due to the fact that in the 
inelastic threshold region (see argumentation above) where the change of the delta potential 
has the biggest effect, the multi-step and elastic scattering dominates the cross section. To 
study the influence of the delta potential we have to compare with experimental data which 
are mainly dominated by the excitation of deltas, i.e. delta mass distribution from vr + 
nucleon invariant mass spectra (see for example Ref. |^) in (p,x)x'7r reactions. 

Concerning the dependence on the angular distribution of the adopted elastic nucleon- 
nucleon cross section we find strong dependence on the adopted parametrisation. Since 
the OSMC calculation is much less time consuming we studied the cross section dependence 
within the OSMC. In Fig. ^ we show the result of the OSMC calculation for isotropic elastic 
nucleon-nucleon cross section (solid line) and for the isospin-independent Cugnon parametri- 
sation [m for the elastic nucleon-nucleon cross section. Since the isospin-independent Cugnon 
parametrisation is strongly forward peaked for all nucleon-nucleon channels, at O = 15° the 
QEP is largely overestimated and at larger angles the QEP is underestimated (Fig. |T2 
(dashed line)). In case of an isotropic nucleon-nucleon cross section the opposite is the case: 
The QEP at G = 15° is much better reproduced, but for all other angles there is a large 
overestimation of the QEP (Fig. |1^ (solid line)). 

If we would do AMD calculations with the isotropic or the isospin-independent parametri- 
sation of the nucleon-nucleon cross section, neither choice would improve the results, which 
we achieved using the isospin-dependent parametrisation of Cugnon. This is a satisfying 
result, since the isospin-dependent parametrisation of Cugnon gives the best reproduction 
of the experimental data among the dicussed parametrisations. 

The strong dependence on the angular distribution of the nucleon-nucleon cross section 
and the effect of multi-step scattering discussed above show how difficult it is to explain the 
considered experimental data for the (p,p') reaction which at first view seem to be very dull, 
uninteresting and easy to understand. 

With the OSMC we checked the dependence of the results on the adopted momentum 
distribution of the target. In Fig. ||, besides the original OSMC calculation (solid line), also 
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calculations with a box momentum distribution of the target (dotted line) are shown. The 
higher momentum components in the target lead to a broadening and reduction of the QEP. 
But in general the (p,p') reaction is not sensitve to the details of the momentum distribution 
in the target. The large variation of width of the QEP between Thomas Fermi and box 
momentum distributions is caused by a large difference of the distributions themselves. 



We did calculations with different mass distributions for the delta discussed in Sec. [II D 
but these only have a minor effect on the result of the calculation. 

We have checked the effect of Pauli blocking in the OSMC. All OSMC calculations 
presented in this paper incorporate Pauli blocking. In the AMD the Pauli blocking is 
incorporated naturally (see Ref. [|r^ for details). We found that only the high momentum 
tail of protons at = 15° is influenced by the consideration of Pauli blocking. 

In general we see in the multi-step decomposition of the AMD calculation in Fig. ^ 
that the contribution of elastic scattering is overestimated. This is especially prominent for 
6 = 30° but also for other angles the defect can be seen. As a first attempt to cure this 
problem we did calculations in which we reduced the elastic nucleon-nucleon cross section 
artifically by 50%. This leads to the results shown in Fig. ^ (dashed line). The QEP can 
be reproduced much better and also the general agreement between experimental data and 
AMD calculation gets better. At B = 40° the smooth momentum dependence of the cross 
section is well described but the cross section is slightly underestimated. 

The reduction of the elastic cross section is also motivated by the fact that we find a 
total reaction cross section of 307 mb for the ^^C(p,p') reaction at i?iab = 800 MeV in the 
AMD calculation, which is too large compared to the experimental results (see Sec. |V^). 

With the OSMC we did calculations with a 70% reduced elastic cross section. For this 
choice of the cross section we have to use Xfactor = 190 mb which is less than the measured 
total cross section. The results are shown in Fig. ^ (dashed line). Also in case of the OSMC 
the reduction of the cross section leads to a better agreement with the experimental data 
for O = 15° and O = 30°. But in the case B = 40° and B = 60° the agreement with the 
experimental data gets worse. 

One possible physical effect which will lead to a reduction of elastic scattering events is 
the correct treatment of the momentum dependence of the mean field. As written above the 
adopted Gogny force does not give the correct momentum dependence of the nucleon-mean 
potential. Instead of [/iv(800MeV) ^ 50 MeV the Gogny force leads to 1/iv(800MeV) ^ 
MeV. Using a repulsive potential would lead to more nucleons diffracted to larger angles. 
This also implies that nucleons would have less interactions with the target since they are 
diffracted from their way through the target. This effect is especially true for nucleons with 
large impact parameter which contribute strongly to the elastic scattering peak. In a future 
investigation this idea should be tested. 

Finally we shortly want to comment on the differences of the AMD and OSMC calcu- 
lations for the ^^C(p,p') reaction at i^iab = 800 MeV from the PWIA calculations by Y. 
Alexander et al. |^6| for the same process. Y. Alexander et al. take one-step elastic and 
one-step inelastic scattering into account. Unlike in the OSMC results discussed in this 
paper the smooth dependence on the outgoing proton momenta for larger angles can be 
better described in their PWIA calculation. This is because they use an unrealistic high 
pion production cross section. Therefore they need no multi-step contributions to get the 
smooth behaviour of the cross section for large angles as we find in the AMD calculations. 
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If one would add the multi-step contributions to their PWIA calculation this would lead to 
a large overestimation of the experimental data. 

Another defect of the PWIA calculations of Y. Alexander et al. is that they fit the heights 
of the QEP in the double differential cross section with an angle dependent scale factor which 
hides the strong dependence on the angular distribution of the elementary nucleon-nucleon 
cross section. 

In the region of phase space where multi-step contributions determine the cross section 
Y. Alexander et al. underestimate the inelastic scattering cross section. This is in good 
agreement with the findings presented in this paper using the OSMC model. 

V. SUMMARY 

In this paper we discussed how to incorporate delta degrees of freedom into the AMD. 
Especially we discussed the inclusion of inelastic channels in the collision term and the 
form of the delta-nucleon potential. For the dclta-nucleon potential we used a Skyrme type 
potential which gives Ua{po) ~ — 30MeV for infinite nuclear matter calculation, thereby 
reproducing the experimental data found in pion-nucleus scattering. 

With the help of the frictional cooling method we calculated delta separation energies 
for 'delta nuclei'. There we found that the minimum energy of the delta-nucleus system in 
the present calculation is achieved when the delta and nucleus wave functions are in the 
state of maximum overlap. Therefore we found that the clustering of AMD wave function 
is refiected in the delta separation energy. 

For a detailed discussion of the (p,p') reaction we developed a one-step Monte Carlo 
(OSMC) model for this reaction. As required by the energy region considered in this paper 
(Eiab ~ 800 MeV) we incorporated both nucleon and delta degrees of freedom using the 
same cross section as done in the AMD. 

As first application of the OSMC we investigated the question of Lorentz invariance for 
nucleon-nucleus scattering for energies E'lab ~ 800 MeV. We found the interesting result 
that, applying a combination of Lorentz transformations, also a Gallilei invariant theory can 
be used to treat nucleon-nucleus scattering in this energy region. Also the calculation of 
in a Gallilei invariant theory can be performed in an approximate Lorentz invariant way. 

Before applying the AMD calculation to the (p,p') reaction, we discussed how to deter- 
mine the width parmeter u of the single particle wave function for outgoing nucleons. The 
determination was made by comparing 'one-step AMD' and OSMC calculations. This is a 
unique way to get parameter independent results in the AMD. 

Both, with extended AMD and OSMC, we preformed calculations of double differential 
cross sections da^ / dp' / dft for the ^^C(p,p') reaction for angle and momentum-dependent 
cross sections. We found a quantitative difference between AMD and OSMC calculations 
which is mainly due to multi-step contributions to the cross section, which are absent in 
one-step calculations. 

Generally we find that OSMC and AMD can reproduce the qualitative behaviour of the 
experimental data. Both the delta peak and the QEP can be seen in the calculation and the 
angle dependence of the cross section is reproduced rather well. In case of the AMD we find 
a quantitative good description of the experimental data expect for the overestimation of the 
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QEP in forward direction. Therefore the agreement with the experimental data improves 
if we perform calculations with a reduced elastic cross section. One reason for the reduced 
elastic cross section could be that the correct treatment of the momentum dependence of 
the nucleon mean field at this energy region leads to a repulsive potential and, because of 
this, to more diffraction of the incoming protons to larger angles. This then leads to less 
elastic scattering probability. This effect would be most serious for big impact parameters. 

For a better understanding of the reaction process we preformed a multi-step decomposi- 
tion of the cross section. To get the smooth momentum dependence seen in the experimental 
data for = 30° and 9 = 40°, a complicated interplay of elastic and inelastic channels is 
important which is approximately reproduced by the AMD calculation. 

One interesting result of the multi-step decomposition of the cross section is that two 
and higher step contributions govern the high momentum tail of cross sections for larger 
angles. The unsatisfying point is that an overall reduction of the elastic scattering cross 
section, which seems to be needed to reproduce the experimental data better, leads to an 
underestimation of multi-step processes and therefore to an underestimation of the cross 
section in the region of phace space where these processes determine the cross section, such 
as the region of the high momentum tail at larger angles. Again the above offered cure to the 
overestimation of the elastic scattering, a correct treatment of the momentum dependence 
of the mean field, would mainly reduce the one-step contributions and not the multi-step 
contributions. 

As one important physical property, which is often discussed but not well understood, 
we discussed the influence of the delta potential on the (p,p') reaction. We found that the 
results for the (p,p') are not very sensitive to the choice of the delta potential. The reason 
for this is that the contributions to the cross section which change due to different delta 
potentials are small and are hidden below multi-step and elastic step contributions to the 
cross section. We pointed out that the dependence on the delta potential will be more 
apparent in calculations for experimental data, where the excitation of the delta dominates 
the reaction, i.e. delta mass distribution from vr -|- nucleon invariant mass spectra (see for 
example Ref. [p5|] ). 

Commenting on prevoius calculations we found that the AMD model gives a much more 
realistic picture of the ^^C(p,p') reaction at i?iab = 800 MeV than the PWIA calculations of 
Y. Alexander et al. in Ref. |^ for the same process. 

Finally we conclude with an outlook: This paper leaves many open questions and ways 
to extend the field of study. First of all a nucleon-nucleon potential should be adopted which 
gives the correct momentum dependence of the nucleon mean field and then the statements 
made about the change of the heights of the QEP should be tested. 

Also the full pion dynamics should be incorporated in the AMD, which is expected to 
have a minor effect on the calculations presented in this paper where we treat a light target 
^^C, since we find that the AA^ — >• NN process is negligble. But for the study of the (p,p') 
reaction for heavier targets the correct treatment of pion absorption (AA^ —>■ NN reaction) 
and pion reabsorption [ttN A) will be much more important. 

With fully incorporated pion dynamics detailed studies of pion production and pion- 
nucleus reactions can be done. 

Another way to follow is to perform a more detailed study of dependence on delta and 
nucleon potential. Especially a momentum-dependent delta potential could be introduced 
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and the difference between momentun-dependent and momentum-independent potentials 
could be discussed. For this investigation we have to compare to experimental data which 
are more sensitive to the adopted potentials. Also we have to introduce a self consistent 
treatment of the delta mass in the AMD for the the NN — > NA process as discussed in Sec. 

Further, (p,n) reactions should be investigated in detail. For this study the the difference 
of transverse and longitudonal delta excition should and can be taken into account in the 
cross section. This will give a nice tool to study the shift of the delta peak in the (p,n) 
reaction. 
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FIGURES 



FIG. 1. Basic Feynman diagram for pion production in nucleon-nucleon collisions. 

FIG. 2. Delta binding energies for various nuclei calculated in AMD with frictional cooling 
method. 

FIG. 3. OSMC of the p(p,n)p7r+ process at E'lab = 830 MeV with different delta mass distri- 
butions. Experimental data are from Ref. (See text for details.) 

FIG. 4. The total reaction probabilities in the OSMC simulation which are determined by the 
total cross sections for proton incident energy E'lab = SOOMeV. 

FIG. 5. One-step scattering simulation for the quasi elastic process for @ = 15° and O = 60° 
at i?iab = 800 MeV with full Lorentz invariant kinematics (solid line) and combination of Gallilei 
and Lorentz invariant kinematics (OSMC) (dashed line). (See text for details.) 

FIG. 6. One-step scattering simulation for the quasi elastic process for = 15° and O = 60° at 
-E'lab = 800 MeV with full Lorentz invariant kinematics (solid line) and Gallilei invariant kinematics 
starting with correct energy (dashed dotted line) and correct momentum (dashed line). (See text 
for details.) 



FIG. 7. The distribution of ^/s in one-step scattering calculations at i?iab = 800 MeV with 
full Lorentz invariant kinematics (solid line) and combination of Gallilei and Lorentz invariant 
kinematics (OSMC; dashed line). (See text for details.) 



FIG. 8. OSMC calculations described in Sec. Ill in comparison with experimental data for 
d?a/dQ,/dp' in ^^C(p,p') reactions for various angles = 15°, 30°, 40° and 60° in the laboratory 
system and £^iab = 800 MeV. Calculations are for free elastic nucleon-nucleon cross section (solid 
line), 70% reduced elastic nucleon-nucleon cross section (dashed line) and box momentum distri- 
bution for the target (dotted line). For the solid and dashed lines the Thomas Fermi momentum 
distribution for the target is used, and for the dotted line free elastic nucleon-nucleon cross section 
is used. (See text for details.) 



FIG. 9. AMD calculation with free elastic nucleon-nucleon cross section (solid line) and 50% 
reduced elastic nucleon-nucleon cross section (dashed line) in comparison with experimental data 
for cfa/dQ,/dp' in ^^C(p, p') reactions at Siab = 800 MeV. Experimental data are the same as in 
Fig. |. 

FIG. 10. Multi-step decomposition of the original AMD calculation (solid line) given in Fig. ^. 
(See text for details). Experimental data are the same as in Fig. ^. 
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FIG. 11. Original AMD calculation (solid line), 'one-step AMD' calculation (dashed line) and 
AMD calculation with shifted delta potential (dotted line) in comparison with experimental data 
for d'^a/dO./dp' in "'^^C(p, p') reactions at -Eiab = 800 MeV. (See text for details). Experimental 
data are the same as in Fig. ^. 

FIG. 12. OSMC calculation with isotropic nucleon-nucleon cross section (solid line), forward 
peaked nucleon-nucleon cross section (dashed line) and box momentum distribution for the target 
(dotted line) in comparison with experimental data for d'^a/dO./dp' in ^^C(p, p') reactions at 
E\ah = 800 MeV. (See text for details). Experimental data are the same as in Fig. ^. 
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